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FLUID FLOW DYNAMICS UNDER LOCATION UNCERTAINTY 

ETIENNE MEMIN* 

Abstract. 

We present a derivation of a stochastic model of Navier Stokes equations that relies on a decomposition of the 

velocity fields into a differentiable drift component and a time uncorrelated uncertainty random term. This type 

of decomposition is reminiscent in spirit to the classical Reynolds decomposition. However, the random velocity 

fluctuations considered here are not differentiable with respect to time, and they must be handled through stochastic 

calculus. The dynamics associated with the differentiable drift component is derived from a stochastic version of 

?-H , the Reynolds transport theorem. It includes in its general form an uncertainty dependent "subgrid" bulk formula 

t-^H that cannot be immediately related to the usual Boussinesq eddy viscosity assumption constructed from thermal 

■^^ molecular agitation analogy. This formulation, emerging from uncertainties on the fluid parcels location, explains 

with another viewpoint some subgrid eddy diffusion models currently used in computational fluid dynamics or in 

0^ ' geophysical sciences and paves the way for new large-scales flow modelling. We finally describe an applications of 

our formalism to the derivation of stochastic versions of the Shallow water equations or to the definition of reduced 

order dynamical systems. 

£>V 1. Introduction. For several years there has been a deep and growing interest in defin- 

ing stochastic models for climate or geophysical sciences J25l |39l (see also the thematic 
^J | issue 13TI . and references therein). This interest is strongly motivated by the necessity of 

devising large-scale evolution models that allow taking into account statistical descriptions 
rJ. J of processes that cannot be accurately specified - to keep, for instance, an affordable compu- 

tational time. This includes small-scale physical forcing (cloud convection, boundary layer 
*ry} ' turbulence, radiative interaction, etc.) and numerical hypothesis and choices operated at the 

£>-». dynamics modeling level, such as a scale coarsening in a given direction, or the selection 

1-jh ' of a particular functional space to express the solution. The modeling simplifications and 

the different unresolved processes involved introduce de facto errors and uncertainties on the 
constitutive equations of the state variables, which are, in general, so complex that only a 
probabilistic modeling can be envisaged. 

In traditional large-scale modeling, the interaction between the resolved flow compo- 
^•s nent and the unresolved components lies essentially in the constitution of a so-called subgrid 

£f~) \ stress tensor, which is usually not related to uncertainty or to an error concept. We believe 

l/~) ■ that it is important to extend this notion in order to take into consideration in a more appropri- 

tH- J ate way the action, on the resolved component, of stochastic processes modeling errors and 

C^ ■ uncertainties of the state variables' dynamics. 

The modeling and the handling along time of such uncertainties are crucial, for instance, 
for ensemble forecasting issues in meteorology and oceanography, where ensembles of runs 
are generated through randomization of the dynamics' parameters, accompanied eventually 
by stochastic forcing that mimics the effect of unresolved physical processes J39l . Problems 
related to the underestimation of the ensemble spread and a lack of representativity of the 
subgrid stress in terms of the model errors constitute in particular problematic limitations of 
the ensemble techniques' forecasting skill. Such stochastic evolution models are also needed 
for data assimilation procedures defined through Monte Carlo implementation of stochastic 
filters referred to as ensemble filters or particle filters in the literature |9][T7). In all these 
situations, a large-scale flow evolution description that includes stochastic forcing terms and 
an uncertainty related subgrid stress expression hence appears necessary. The root problem 
of this general issue lies essentially in the construction of adequate stochastic versions of the 
Navier-Stokes equations. 
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2 E. Memin 

The establishment of sound stochastic dynamical models to describe the fluid flow but 
also the evolution of the different random terms encoding uncertainties or errors is a difficult 
issue compounded with an involved mathematical analysis. As initiated by |2| and intensively 
studied by many authors since then (see for instance the review ifPZl . and references therein) 
such a formalization has been mainly considered through the addition of random forcing 
terms to a standard expression of Navier-Stokes equations. However, this construction is 
limited by a priori assumptions about the noise structure. Furthermore, the question whether 
this noise should be multiplicative or additive immediately arises. 

In this work, we wish to explore the problem via a reverse strategy. Instead of consider- 
ing a given - eventually simplified - dynamics and then to supplement it with random forcing 
terms to model errors carried by unresolved or unknown processes, we will start from a gen- 
eral Lagrangian stochastic description that incorporates uncertainties on the fluid motion. The 
sought-after Eulerian dynamics is then deduced from this general stochastic velocity descrip- 
tion and standard physical principles or approximations. This construction, reminiscent to 
the framework proposed by lETl . has the great advantage to let naturally emerge determinis- 
tic and stochastic uncertainty terms related to the different errors transported by the evolution 
model. Deterministic approximations or stochastic implementations of the dynamics can then 
be considered on solid grounds. 

Following this route, we aim in this paper at devising stochastic dynamics for the descrip- 
tion of fluid flows under uncertainties. Such uncertainties, modeled through the introduction 
of random terms, allow taking into account approximations or truncation effects performed 
during the dynamics constitution steps. This includes for instance the modeling of the unre- 
solved scales interaction in Large-Eddies Simulations (LES) or Reynolds Average Numerical 
Simulations (RANS). These uncertainties may also encode numerical errors, unknown forc- 
ing terms or uncertainties on the initial conditions. They gather the effects of processes one 
does not wish to accurately model. However, as they propagate along time and interact with 
the resolved components their effects must be properly taken into account. 

We will assume throughout this study that the fluid particles displacements can be sep- 
arated in two components: a smooth differentiable function of time and space and an un- 
certainty function uncorrected in time but correlated in space. This latter component is for- 
mulated as a function of Brownian motion and the whole displacement is defined as an Ito 
diffusion of the form: 

dX(x,t) = w(X(x,t),t)dt + cr(X(x,t),t)dB t , (1.1) 

where X : il x R + — > O is the fluid flow map, which represents the trajectory followed 
by a fluid particle starting at point Xu—o(x) — x of the bounded domain 0. This consti- 
tutes a Lagrangian representation of the fluid dynamics and dX(x, t) figures the Lagrangian 
displacement map of the flow at time t. In expression ( II. lb . the velocity vector field, w, 
corresponds to the smooth resolved velocity component of the flow. It is assumed to be a 
deterministic function (of random arguments) and to have twice differentiable components: 
(u, v, w) £ C 2 ($7,K). When it does not depend on random arguments this velocity field 
represents the expectation of the whole random velocity field. This is the situation encoun- 
tered in the case of mean field dynamics. The second term is a generalized random field that 
assembles the unresolved flow component and all the uncertainties we have on the flow. The 
combination of both velocity fields provides an Eulerian description of the complete velocity 
fields driving the particles: 

U(x,t) =w(x,t)dt + cr(x,t)dB t . (1.2) 

The whole random field, U(x, t), should be a solution of the Navier-Stokes equations and 
is defined here as the combination of a stochastic uncertainty component and a deterministic 
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"resolved" component driven respectively by unknown characteristics, er and w that have to 
be determined or specified. 

This framework is also related to the work of [j5]. Nevertheless, their study aimed at 
building a Monte-Carlo Lagrangian representation of the Navier-Stokes equations. It is lim- 
ited to a constant diffusion tensor a and exhibits additional difficulties to deal with boundary 
conditions \1\. In this work, we do not look for a stochastic Lagrangian representation of the 
deterministic Navier-Stokes equation. We seek instead, a deterministic smooth representa- 
tion of the drift velocity component corresponding to a solution of a stochastic Navier-Stokes 
equation. The goals are thus in some way opposite. Our objective can be interpreted in terms 
of Large Eddies Simulation (LES) or Reynolds Average Numerical Simulation (RANS), as 
we aim at separating a resolved flow component from an unresolved one, respectively de- 
fined as a differentiable component and a time uncorrelated uncertainty component. From 
this prospect, the approach proposed here is closer in spirit to RANS techniques than to the 
LES paradigm as no spatial filtering is applied (see j22l [26l [36), for extensive reviews on 
the subject). Furthermore, let us point out that the uncertainty component lives at all the 
hydrodynamical scales. Thus, there is no spatial scale separation principle here but rather a 
temporal decomposition in terms of a highly oscillating process with no time differentiation 
property and a smooth differentiable component. This appraoch is also close in spirit to the 
separation in term of a "coherent" component plus noise operated through adaptive wavelet 
basis iflOlfTTI . However, contrary to this approach relying on a Galerkin projection with an 
adaptive scale thresholding, our decomposition makes appear a diffusion tensor assembling 
the action of the unresolved uncertainty component on the resolved component. 

More precisely, the resolved component of the Navier-Stokes equation we consider in- 
cludes in its general form an anisotropic diffusion term that emerges due to the presence of the 
random uncertainty term. By analogy with the Reynolds decomposition and LES, we refer to 
this tensor as the subgrid stress tensor. However, it is important to outline that its construc- 
tion differs completely from the Reynold stress definition. As a matter of fact, in our case 
the unresolved fluctuating component is a non differentiable random process, and stochastic 
calculus differentiation rules have to be used. We will see that in the general case the re- 
sulting subgrid stress cannot be immediately related to the usual eddy viscosity assumption, 
formulated in the nineteen century by Q from thermal molecular agitation analogy - also 
commonly referred to as Boussinesq assumption in the litterature - and that is still intensively 
used in the Large Eddies Simulation paradigm since the work of BUI and l23l . 

The simpler form of eddy diffusion will be recovered only when the uncertainty will be 
confined to homogeneous random fields such as the Kraichnan random field Ifl9ll . For partic- 
ular forms of the uncertainty random component, this formulation will explain with another 
viewpoint some subgrid eddy diffusion models proposed in computational fluid dynamics or 
in geophysical sciences. Following this route, an appropriate definition of uncertainty models 
adapted to given situations allows opening the way for new large scales flow modeling. 

In a very similar way as in the deterministic case, the representation we propose will be 
determined from a stochastic version of the Reynolds transport theorem relying on a specific 
model of the uncertainty random field. 

2. Construction of the uncertainty random field. Before presenting in detail the deriva- 
tion we consider to built stochastic fluid flow evolution models, we need to specify the ran- 
dom component that will encode the uncertainty associated to the fluid parcels location. This 
random field we refer hereafter to as Brownian avatar (as a smooth idealization of a dense 
Brownian map) is defined in the following section. It is built from a finite dimensional dis- 
crete set of standard Brownian variable and tends in law to a continuous limit, which formally 
enables simple calculus. 
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2.1. Brownian motion field avatar. This random field, denoted B t : Q — ¥ R d , is built 
from a finite dimensional discrete set of standard Brownian variables as 

1 " 
B?(a;) = - 7 =VB t (x < )^(x-aJi), (2.1) 

where B t = {Bt(aij), i = 1, . . . , n} is a set of independent d-dimensional (with d = 2 
or 3) standard Brownian motions centered on the points of a discrete grid S = {xi,i = 

1 n}Cfi and <p is a Gaussian function of standard deviation v. It is immediate to check 

that B™ is a zero mean Gaussian process, with uncorrelated increments and whose spatial 
covariance tensor tends for an infinite number of grid points to a limiting covariance: 

Q = lim E[B?(x)Br(y)] = ttp^ v (x - y)l d . (2.2) 

As B™ is a Gaussian process, it tends in law to a zero mean continuous process with the 
same limiting covariance. This limiting process will be denoted in a formal way through a 
convolution product: 

Bfktp v (x)= / B t (x')ip u (x - x')dx' , (2.3) 

and will be identified to B t in the following. The covariance operator Q is positive definite 
and of finite trace 



\B\\l = ^ 



V/ ipl(x-Xi)dx = td{ATtv 2 )- d/2 . (2.4) 

, JR d 



Hence, the energy of the Brownian avatar does not depend on the number of grid points used 
for its construction but only on the standard deviation of the Gaussian smoothing function. 
At finite time horizon, when the standard deviation tends to infinity, the energy brought by 
the uncertainty tends thus to zero and the SDE (II . lb behaves almost like a deterministic 
evolution law. At the opposite, in the zero limit of the smoothing function standard deviation, 
the variance - or the energy - becomes unbounded and the Ito stochastic integral can be no 
more defined. For non zero standard deviation the covariance operator is of finite trace and 
the Brownian avatar constitutes a Q- Wiener process 11331 . 

2.2. White noise avatar. The analogue of the white noise on the bounded domain H is 
similarly defined in a formal way from the generalized function dB t = dB t * tp v and, a t , a 
linear bounded deterministic symmetric operator of the space of mean square integrable func- 
tions with null condition outside the domain interior K. d /f2 U d£l. The uncertainty component 
reads then 

<r(x,t)dB t = I <T t (x,y)dB t (y)dy. (2.5) 

This operator <x t is referred in the following to as the diffusion tensor. The covariance tensor 
of the uncertainty component <j t dB t , (also called the two points correlation tensor) reads: 

1 " 

Q = lim —dt6(t — s)'S^<T(x t ;t)*ip v (xi)<r(;y,t)*(pv(xi) 

n— »oo n *• — ' 

i=0 

= dt5(t- s)cr Vv {x,t)cT v>v (y 1 t), (2.6) 



Fluid flow dynamics under location uncertainty 5 

where o~ v (x, y, t) — <r(x,*,t) * ip v (y) is a filtered version of cr along its first or second 
component as it is symmetric. Processes of central importance in stochastic calculus are the 
quadratic variation and covariation processes. They correspond to the variance and covari- 
ance of the stochastic increments along time. The quadratic covariation process denoted as 
(X, Y) t , (respectively the quadratic variation for Y = X) is defined as the limit in proba- 
bility over a partition {ti, . . . , t n } of [0, t] with t\ < ti < ■ ■ ■ < t n , and a partition spacing 
5ti = U — ti-i, noted as \5t\ n = m&xdti and such that |<5tL — ^ when n — > oo: 

i 
p n—l 

(X,Y) t = lim T t {X(t i+1 )-X(U))(y(U ¥1 )-Y{t i )) T . 

' ' n_ *' i=0 

For Brownian motion these covariations can be easily computed and are given by the follow- 
ing rules (B, B) t = t, (B, h) t = (h, B) t = (h, h) t = 0, where h is a deterministic function 
and B a scalar Brownian motion. In our case, the quadratic variation of the uncertainty com- 
ponent reads: 

ft rt 



(f (*(x,t)dB t y, f (*(x,t)dB t 
Jo Jo 

lim — / } } <x lfe (a:,», s) * (pv(xg)cr k; > (u,x,s) *ip v (xAds 



/0 1=0 k 



I 

V)cr*(a! J «)o*{(»,*)dfl 

k 
t 



I a ll {x,s)ds a.s. (2.7) 



A 

Its time derivative corresponds to the diagonal of the spatial covariance tensor. We can note 
that for homogeneous random fields cr(x — y) these expressions simplify. As a matter of fact, 
in this case, the random component cr(x, t)dB t = <r(», t) * dB t (x) has a covariance tensor 
that is defined as 

Q = dtS(t - s)(t(; t) * <r(«, t) * (fi^^x - y). (2.8) 

It remains thus homogeneous and its quadratic variation tends to a spatially constant d x d 
matrix for an infinite number of grid points: 

1 ™ f 

lim — \ \cr(*,t) * (p„(x — Xi)] 2 dt =dt / \(r(»,t)-k(p v (x)] dx 

n-^oo n Z-j J n 



i=0 



=T,{t)dt a.s. (2.9) 



Besides, we remark that spatial derivatives of the white noise avatar can be defined and ma- 
nipulated in a simple manner: 

1 - 

d Xt {cr(x)dB t ) = lim — = S^ d Xi cr ip ^(x,Xj)dB t (xj) 

V 3=0 

=d Xi cr Vu {x)dB t . 

As a direct consequence, we note that a divergence free random component is thus necessary 
associated to a divergence free diffusion tensor. The derivative operator of the identity diffu- 
sion tensor cr(x, y) = S(x — y) is also still well defined for a nonzero standard deviation and 
acts solely on the white noise avatar Gaussian smoothing function. 
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2.3. Kraichnan turbulent model. Toys models pioneered by [|T9l and intensively ex- 
plored for the theoretical analysis of passive scalar turbulence |[T4l [191 l24l can be easily 
specified with such a model. The Kraichnan random field model 

dt c t =V*i>Z*f C *dB t , (2.10) 

is formally defined from the divergence free projector V, a power function f^(x) = C<j 1 1 a; 1 1 ^/ 2 
for < C < 2 and a band-pass cut-off function, ipj, on the so-called inertial range that lies 
at high Reynolds number between the short dissipative scale l D = 1/k and the large inte- 
gral scale L = I/7 at which the forcing takes place. The constant Ct has a dimension of 
£ 1 ~z/t, which follows from dimensional analysis (the energy transfer rate e = Ef^) ~ |y 
and e ~ C? ^-)- The spectral correlation tensor of dQ is defined as: 

Qikh = m-t-^ - 1^)(^) 2 , 

and the variance (or time derivative of the quadratic variation process) of this homogeneous 
random field reads 

d{${x),$(x))=dt [[V*K*f<(x)] 2 dx a.s, (2.11) 



which is constant and given for a band-pass filter, ipj, defined as a box filter as 

dtC c d - 1 2ir d/2 
(2w) d d T(|) 



tt%nrCHL< -*„)*. (2-12) 



The square root of this quantity represents the mean absolute value of the velocity field com- 
ponents. It blows up when the integral scale tends to infinity and exhibits a strong depen- 
dency on the largest length scale. This injection scale dominates the turbulent energy and 
corresponds to the distance beyond which velocities decorrelate. 

We will see that such random fields, which are both homogeneous in space and uncor- 
rected in time, lead to a simple form of the dynamics drift term involving an intuitive eddy 
viscosity term defined from the uncertainty random field variance. However more general 
random fields will let rise a formulation of the Navier-Stokes equations in which emerges an 
anisotropic diffusion term that cannot be immediately related to the usual eddy viscosity as- 
sumption first formulated in the nineteen century by [5 L| and that is still predominently used 
in the Large Eddies Simulation paradigm. 

3. Stochastic Reynolds transport theorem. In a similar way as in the deterministic 
case, our derivation relies essentially on a stochastic version of the Reynolds transport theo- 
rem. The rate of change of a scalar function q within a material fluid volume V(t) transported 
by dl.ll i is: 

d fq(x, t)dx = [{dq t + [V-(qw) + h\V-*\\ 2 q - ^ \^-^~ (Otfff)|^ =0 ]dt 
JV(t) JV(t) Z ~J l OXiOXj 

+ V-(qcrdB t )}dx. (3.1) 

In this expression the third term must be computed considering the diffusion tensor is diver- 
gence free, and the tensor a(x) denotes the uncertainty variance defined in (12.7) . 



'referred sometimes for that reason as the Boussinesq assumption. 
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To give some insights on how this expression is obtained, let us consider a C 2 (ft) func- 
tion g : SI — >• K that vanishes outside V(t) and on dV(t). Let us also define a scalar function 
</> on fi x K + such that: 

<f)(X t ,t) = g(x ). 

Let us outline that this function cannot be a deterministic function. As a matter of fact, if it 
was the case, its differential would be given by a standard Ito formulae: 

dfat, X t ) = {d t (j> + V(f> T w + \Y1 d ( X l X i)-ff-^-) dt + V(f> T <rdB t . 

This quantity cancels as </> is transported by the flow. A separation between the slow deter- 
ministic terms and the fast Brownian term would yield a very specific uncertainty liying on 
the function iso-values surfaces ( V(/) T er = 0), which is not the general goal followed here. 

For a fixed grid point, y, 4>(Xt , t) is hence necessarily solution of a stochastic differential 
equation of the form: 

dfa{y) = v(y, t)dt + f(y, t) T dB t . (3.2) 

As fax, t) is a semi-martingale, the differential of <fr(X tl t) involves the composition of two 
stochastic processes. Its evaluation requires the use of a generalized Ito formulae usually 
referred in the literature to as the Ito-Wentzell formulae (theorem 3.3.1 JpTI). This extended 
Ito formulae incorporates in the same way as the classical Ito formulae^ quadratic variation 
terms related to process Xt but also covariation terms between Xt and the gradient of the 
random function fa. Its expression is in our case given by: 

dfat, x t ) =dfa + wdXt + J2 d (^r> (<rB t y)dt + ±j2 d ( x i x i) -^t^ dt 

i i,j ? 3 

=0. (3.3) 

For a deterministic function or for uncertainties directed along the surface of function iso- 
values it can be immediately checked that the standard Ito formula is recovered. The quadratic 
variation term involved in (13.3) is given through (12.7) as 

d(Xt,Xi)=a ij (x,t), (3.4) 

and the covariation term reads: 

d{—± 1 (crBt) l )= lim - V V(- 5 — (x,»,t) * <p v {xe))<T %: > (x,;t) *(p v {x t ), 

OXj n->oo n * — ' * — ' OX; 



From these expressions and equations (13.2H3.3) we infer: 

f(x,y 7 t) T = -Vfax,t) T <T(x,y,t), 

d 2 d> _ in ,dtr'P 



v(x,t) = -V(^ T w + J2^a 



Vo 



'-,/ 



2 3 dxidxj ' dxi 



H-*n 



2 relevant only to express the differential of a deterministic function of a stochastic process 
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and finally get: 






dffo = Ccj>dt-Vcj> T (TdB tl 




1,3 



(3.6) 



This differential at a fixed point, x, defines the equivalent in the deterministic case of the 
material derivative of a function transported by the flow. The differential of the integral over 
a material volume of the product q<fi is given by 



/ q<j>(Xt,i)dx — d / q<j>dx, 
Jv(t) Jn 

(dq t (j> + qd(j)t + d(q, <p)tdx, 



where the first line comes from <fi(t,x) = if a; £ Cl\V(t) and the second one from the 
integration by part formula of two Ito processes. From (13. 6) this differential is hence: 

/ (dq t 4> + qC4> + V 4>aV q)dtdx - I qV<\> T crdB t . 
Jn Jn 

Introducing the (formal) adjoint C* of the operator C, this expression can be written as 

(dqt + C*q-V-(aVq))dt+V- (qcrdB t )) §dx. 

With the complete expression of C* and if <f>(x,t) — > Uv(*)/9V(t)> where II stands for the 
characteristic function, we get the extended form of this differential: 



f[dq t + {V-(qw)-J2^-(q< 3 ^) 

JVM ;; OXi 



/V 0) ^ 



\ E dx^x~ {aijq) " E ^a^g - ]T a*dl iX .q}dt + V-{qadB t )]dx, 



IJ ■> I] 

which can be further simplified as 

P 
z Ijk vXibx 



f[dq t + V-(qw)dt - \ £(a*cr^-^- + d Xj a?d Xi a k Jq+ 

JV{t) A ljk OXiOXj 

2<jfd Xi a h Jd Xi q- d Xl ul k d X] a^q)dt + V-(q<rdB t ))dx. 

At last this expression may be more compactly written as: 

f Id 2 1 

J(dq t + V-{qw)dt - 2 E dx-dx- ( aij9 ^'=° dt + 2 11 V ' <T ^ qdt + V-(q<rdB t ))dx, 



where the third term must be computed by considering a divergence free diffusion tensor. 
Let us note it is now straightforward to get the expression of the transport theorem for a 
divergence free turbulent component: 

d q(x,t)dx = / dq t + [Vqw - -^ df ' j (a i jq)]dt + VqcrdB t dx. (37) 

Jv(t) Jv(t) l*r< 
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4. Mass conservation. This expression of volumetric rate of change relation allows us 
stating the mass conservation principle under fluid particles location uncertainty. Applying 
the previous transport theorem to the fluid density p(x, t), we get a general expression of the 
mass variation : 

/l d 2 1 

[dpt + (Vipw) - - Y, a^r^Wo + 2 II V - (T W 2 P) dt + V-{ P adB t ))]dx. 



(4.1) 



A mass conservation constraint on the transported volume implies cancelling this expression. 
Besides, as the volume is arbitrary this provides the following constraint 

dp t + V-(pw)dt = i(£ ——( aij p) w<r ^ - ^\\V-a\\ 2 p)dt - V-{padB t ). (4.2) 

4.1. Incompressible fluids. For an incompressible fluid with constant density, can- 
celling separately the slow deterministic terms and the highly oscillating stochastic terms 
of the mass conservation constraint (14.2), we obtain: 



V-(<rdB t )=0, (4.3) 

1 y^ d 2 a,ij 

2 4r* dxidxj 

Considering we are seeking a volume preserving flow for w, this system boils down to an 
intuitive system of incompressibility constraints: 

V-(crdB t ) = 0,V-w = 0, (4.5) 

coupled with a non intuitive additional constraint on the quadratic variation tensor: 

V-(V-a) = 0. (4.6) 

For the Kraichnan model (or for any divergence-free homogeneous random fields) this last 
constraint is naturally satisfied, as its quadratic variation is constant and the system comes 
down to the classical incompressibility constraint. In the same way, any homogeneous divergent- 
free random fields associated to a volume preserving drift leads to mass conservation. 

4.2. Isochoric flows and isoneutral uncertainty. For divergence-free volume preserv- 
ing flows (refered as isochoric flow) with varying density we get a mass conservation con- 
straint of the form: 

1 d 2 
dpt + Vpwdt - - ^ (pa lj )dt = V padB t . (4.7) 

Zi (J JL %J-'j 

An interesting property emerges if the uncertainty has a much larger scale along the den- 
sity tangent plane than in the density gradient direction. This situation occurs in particular in 
oceanography or in meteorology where the fluid is stratified and the horizontal scale much 
larger than the vertical scale. In a such case, it is essential that the large-scale numerical 
simulations preserve the natural fluid stratification and consequently to define subgrid mod- 
els with controlled diffusion along the iso-density surfaces. This behaviour can be easily set 
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up by defining the diffusion tensor as a diagonal projection operator so that the uncertainty 
crdB t lies on the isodensity surfaces: 

ii xij d Xi p(x)d x p(y) 

(j lJ = 5 13 — — j| 5 (x - y). (4.8) 



Together with the small slope assumption [151 (\/{dxP) 2 +{d y p) 2 « d z p) the diffusion 
tensor and the quadratic variation can then be written as a matrix function 

/ 1 a x (x)\ 

a(x) = 1 a y {x) , (4.9) 

\a x (x) a y (x) \oc(x)\ 2 J 

where we introduced the neutral slope vector a. = (a x , a y ,0) = — (d x p/d z p, d x p/d z p, 0). 
Owing to the divergence free condition of the diffusion tensor, this slope vector is necessarily 
constant along the depth axis d z a x = d z a y = and divergence free V- a. = 0. This yields 
for the mass preserving equation a diffusion along the density tangent plane. The density 
evolves as the deterministic model: 

^ + Vpw = I J2 d Xi ( aij d Xj p), (4.10) 

ij 

where we considered the divergence free constraint V • a = 0. This type of anisotropic 
diffusion corresponds exactly to the so-called isoneutral diffusion currently used to model 
unresolved mesoscale eddies in coarse scale resolution of ocean dynamics simulations 03). 
As mentioned previously, the divergence free uncertainty induces a slope vector that is also 
divergence free and independent of the depth. This enforces further the density surfaces to be 
quasi-harmonic (for smooth depth density variation) in the horizontal planes and provides a 
strong stratification of the density organized as a pile of pancakes with no maxima excepted 
on the domain horizontal boundaries. To our knowledge those additional constraints are not 
taken into account in the isoneutral diffusion and according to our interpretation this comes 
to consider uncertainties that are not volume preserving. 

In the case of the Kraichnan model the density variation is simplified as 

dp t + V pwdt - 7- V] Apdt = VpadB t , (4.11) 

and when the drift term does not depend on random argument (i.e. if it corresponds to the 
flow expectation w = KdX t as in the case of a mean field dynamics) the density expectation 
evolution is a classical intuitive advection diffusion equation: 

dtp + Vpw = -■yAp. (4.12) 

The mass conservation constraint and the stochastic version of the Reynolds theorem 
allows us now expressing the linear momentum conservation equations. 

5. Conservation of momentum. Newton's second law states that the net force acting on 
the fluid is equal to the rate of change of the linear momentum. In our case, the flow evolution 
is described through an Ito diffusion (ll.lt . where the drift component w(x,t) represents 
the unknown flow we wish to estimate whereas the noise term denotes the fluctuating part 
either caused by physical processes that are neglected or generated by numerical errors and 
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coarsening processes. Whatever the type of uncertainties considered on the fluid particles 
location, they introduce de facto an uncertainty on the forces. In order to take into account all 
the uncertainty sources within the momentum equations we consider a stochastic conservation 
principle of the form: 

d p(w(x,t)dt + (r(x,t)dB t )dx = F(x,t)dx. (5.1) 

Jv(t) Jv(t) 

In this equation the acceleration is highly irregular and has to be interpreted in the sense of 
distribution. For every h £ C^°(M+): 

h(t) F(x,t)dxdt = - h'(t) o-(x,t)dB t dxdt+ / h(t) d pw(t,x)dxdt. (5.2) 
Jv(t) J Jv(t) J Jv(t) 

Since both sides of this equation must have the same structure, the forces can be written as: 
h(t) F(x,t)dxdt = - h'(t) a(t,x)dB t dxdt 

JV(t) J JV(t) 

+ h(t) (f(t,x)dxdt + 0(t,x)dB t )dx. (5.3) 
J Jv(t) 

The first terms of the right-hand sides of equations (15.2) and ( 15.3b are identical and cancel out. 
The second right-hand term of equation (15.2) corresponds to the derivative of the momentum 
associated to the resolved velocity component, whereas the second right-hand term of (15.3) 
provides us the structure of the forces under localization uncertainty. We now further develop 
those different terms. 

According to the previous section results, the transport equation applied to the linear 
momentum gives for each component of the velocity: 

d puiidx = l{d{pw i ) t + V-{pw l w) + \\V-(T\\ 2 pw i -^ -- — - — (a jk pwi)i v<T=0 )di+ 
J V (t) Jv(t) ~t 2 9x 3 dxk 

V-(pwiO-dB t ))dx. (5.4) 

As for the forces, we will consider that only body forces and surface forces are involved 
(i.e. there is no external forces). A direct extension of the deterministic case provides us the 
surface forces as 

£= / -V(pdt + dp t ) + p{AU + lv(V-U)), 
Jv "J 

where p is the dynamic viscosity, p(x, t) denotes the deterministic contribution of the pres- 
sure and dp t is a zero mean stochastic process describing the pressure fluctuations due to the 
random velocity component. The gravity force is deterministic and standard. 

Expressing the balance between those forces and the acceleration ( 15.41 ). incorporating 
then the mass preservation principle (14.2) . and finally equating, in the one hand, the slow 
temporal bounded variation terms and, in the other hand, the Brownian terms, we get the 
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expression of the flow dynamics: 



(— +wV w)p- -^a l]P — — -^ — — | w =o 



9tu 



dt 2 ^ dxidx; ^ dx-j dxi 

1,3 J 1,3 

pg-V P + p(Aw + -V{V-w)) 7 (5.5) 
Vdpt =-wV T p ( rdB t +p{A{ ( TdB t ) + ^V{V-{(TdB t ))), (5.6) 



1 d 2 1 

d ft + (V-(pto)--2J— ^— (a l3 p) l ^ =0 + -\\V-<T\\ 2 p)dt = V-{p ( TdB t ). 



dxidxf ljr '^=° 2 1 



(5.7) 



This system provides a general form of the Navier Stokes equations under location uncer- 
tainty. Similarly to the Reynolds decomposition, the dynamics associated to the drift compo- 
nent includes an additional stress term that depends on this resolved component and on the 
uncertainty diffusion tensor. In order to specify further the different terms involved in this 
general model, we will examine simpler particular instances of this system. In the following, 
we will confine ourselves to the case of an incompressible fluid of constant density. Consid- 
ering as a first example the divergence free isotropic Kraichnan model defined in ( 12.10b and 
its associated transport equation, we obtain straightforwardly the following Navier-Stokes 
equation: 

,dw __„ 1 . . _ . 

( — — h w V w — -f-Aw)p = pg — \7p + pAw 

Vdp t = -p{wV T )d£ t + ftAdtt, 

V-to = 0. (5.8) 



The first equation of this model corresponds to a traditionnal turbulent diffusion modeling 
relying on the Boussinesq assumption with a constant eddy viscosity coefficient. Let us note 
that for this system the second equation related to the random pressure term is not needed 
to compute the resolved drift w. The uncertainty is in that case specified a priori as an 
homogeneous random field. Considering a more general divergence free random component, 
we obtain: 



(— +wV T w--^2 d Xl d Xj (a v w))p = pg-Vp + pAw 

Vdpt = -p(wV T )crdBt + pAcrdB t , 

V-io = 0,V-er = 0,V-(V-a) = 0. (5.9) 



This model of Navier-Stokes equations includes now a more general diffusion term that can- 
not be directly related to the traditional Boussinesq eddy viscosity formulation anymore. 
Furthermore, the incompressibility condition on the variance tensor provides a nonlocal con- 
straint on the subgrid model, which is lacking in usual eddy viscosity models[20|. One can 
point out that for divergence-free incompressible uncertainty models this term is globally 
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dissipative. As a matter of fact the energy associated to the subgrid term reads 

-)2 



'W k W k 



f d f d d 

Jn ff dxidxj Jn.dxi ^f d x 

- [ V VwlaVw k - J (V • a) V Vt 

Jn k Jn k 

= [ \\Vw\\l- f(V-a)V(hw\\ 2 ). (5.10) 

Jn Jn L 

The first term || Vio||^ = J^ • a,ijd Xi w-d Xj w, is always non-negative as, a = <Tcr T , is semi- 
positive definite. The second term associated to the incompressibility constraint, V • Vo = 0, 
cancels: 

f (V • a)V(i|M| a )ds = \ f (V • (V • a \\wf))dx (5.11) 

- / ||ii;|| 2 (V ■ a) • nds 
2 Jen 



2 Jas 
0. 



The energy of the subgrid term reduces to 



j WT Y.^Q^.^ W ) dX = - j H V ™!la^ (5-12) 



and is thus globally dissipative. Such model paves the way to set up subgrid models and 
stochastic expression of Navier-Stokes equations for different forms of the diffusion tensor. 
It can be noted that this subgrid model provides the same dissipation as 

f w T y2d Xi {a ij d Xj w)dx = - [ \S7w\ 2 a dx, (5.13) 

Jn , 3 Jn 

which involves a diffusion operator Vf = ^ ■ d Xi (a 1 - 7 ^ /) similar to diffusions introduced 
in LES lfl6l through the concept of spectral vanishing viscosity operator initially introduced 
for ID problems BT1 . No clear consensus, nevertheless, exists in the 3D case on the form 
such an operator should take 



5.1. Link to the Smagorinsky subgrid model. The Smagorinsky subgrid model is 
widely used in large eddies simulation approaches. This eddy viscosity model, built from 
the Boussinesq assumption, imposes that the proper directions of the subgrid stress tensor are 
strictly aligned with those of the resolved rate of strain tensor. The subgrid model derived 
from our uncertainty analysis does not rely on such assumption. Nevertheless, a relation 
with the Smagorinsky model term can be easily emphasized. As a matter of fact, setting the 
quadratic variation tensor to the form: 

a = C\\S\\I, (5.14) 

with C a constant and where ||S|| = |E«(^b 4 u^ + d Xj w' 1 ) 2 ] 1 ^ 2 denotes the Frobenius 
norm of the rate of strain tensor, we get an uncertainty stress tensor term that reads: 

Y,d x AAa^w k ) = Y^d^iWSW&w*) (5.15) 

ij ij 

= 2Y,d Xj \\S\\d X] w k + \\S\\kw k + A||S||«A (5.16) 

3 
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This term complemented by 2 V ■ d Xj WSWdx^ — A||S'||w fc provides the standard trace free 
Smagorinsky subgrid stress, V • r, where 

t = C\\S\\S. (5.17) 

The complementary term may be rewritten as 

2d Xk Y,d Xj (\\S\\)w> -2Y,d X] d Xk (\\S\\)w> - A\\S\\w k . (5.18) 

i i 

We observe that the first right-hand term is a gradient term that can be compensated by a 
modified pressure as it is usually done considering deviatoric expression of the Reynolds 
stress tensor. However, unless assuming that the rate of strain magnitude is a very smooth 
function with flat iso-surfaces there is no particular reason to cancel the two other terms. 
Let us note this condition allows fulfilling the variance incompressibility constraint V • V • 
a = 0. According to our interpretation, the Smagorinsky model constitutes thus an adapted 
uncertainty model only for smooth deformations with flat rate of strain tensor norm. 

In the examples provided previously, the uncertainties have been fixed a priori. In such 
situations the second equation of system ( 15. 9\ is not required unless a realization of the os- 
cillating random pressure field is needed (which is a direct possibility of our formulation). 
However, when the diffusion tensor is not specified the system exhibits another degree of 
freedom and this tensor or at least the quadratic variation tensor must be estimated in order 
to fulfill the balance of Newton's second law. 

5.2. Diffusion and quadratic tensors estimation. The second equation of ( 15. 9t pro- 
vides a mean to proceed to the estimation of the diffusion tensor. Projecting this equation 
on the divergence-free space tensor, for a given increment e = y — x, we obtain a matrix 
screened Poisson equation with variable source functions: 

Aa(x, e + x) = -V * (wV t (t)(x, e + x), (5.19) 

M 

where V indicates the divergence free projector, applied to the column vectors of tt>V T er. 
A symmetric solution can be imposed adding and averaging the transposed problem on both 
sides: 

Acr(x,e + x) = — [V * (wV T a) + (erVw T ) *T](x, e + x). (5.20) 

Formally this equation should be solved for each increment value, which is obviously unreal- 
istic in practice. Nevertheless, with the hypothesis of an homogeneous decorrelation function 
the resolution can be led only for the diagonal (e = 0); the complete tensor cr(x, y) is then 
inferred from those values and the decorrelation function. Another solution consists to as- 
sume a rapid decorrelation and to solve the system for a narrow band around the diagonal. 
An iterative solution can be built assuming an initial diffusion tensor, and solving these equa- 
tions from the current values of w and <x. Let us point out that the quadratic variation tensor, 
a, requires finally an ultimate projection to stick to the last constraint of the system ( 15.9I >. 

6. Shallow water model. We describe in this section an application of our formalism 
to the derivation of a stochastic model of the shallow-water equations. 

The Shallow water equations constitute one of the simplest model that can be used to de- 
scribe the evolution of mean horizontal components of atmospheric winds or oceanic streams. 
This 2D model is valid for problems in which the vertical dynamics can be neglected com- 
pared to the horizontal effects and cannot be used when 3D effects become essential. It is 
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derived from a 3D incompressible model by depth averaging assuming the depth of the verti- 
cal domain is "shallow" compared to the horizontal domain. In order to establish the Shallow 
water approximation, the pressure variable is assumed to follow quite correctly an hydrostatic 
equilibrium relation given by: 

d / = -gp, (6.D 

oz 

where the density is constant on a shallow layer of the fluid. 

In our case, the whole pressure function is given as the summation of a deterministic 
pressure contribution, p', and, dp, a zero mean random pressure: 

P(x, y, z, t) = p'(x, y, z, t) + dp(x, y, z, t). (6.2) 

Assuming the external force is due to gravity and neglecting the friction forces, the hydro- 
static balance for a deterministic system can be regarded as a boundary layer assumption in 
which the vertical acceleration of the fluid is suppose to be null. Considering a system with 
uncertainty under the same hypotheses, that is to say: no friction and a vertical acceleration 
compensated by the subgrid diffusion, we obtain an identical hydrostatic relation (I6.lt . How- 
ever, two different choices are possible. Either the hydrostatic balance can be assume to hold 
for the global pressure or only for the deterministic component. In the first case, in comple- 
ment to an hydrostatic balance on the deterministic pressure component, the random pressure 
component is necessarily such that 

d z d Pt = 0. (6.3) 

A strong hydrostatic assumption on the global pressure component hence requires that the 
random pressure term is constant along the vertical axis. In the second case, there is a sup- 
plementary degree of freedom and only a boundary condition on the random component will 
have to be imposed. 

With the assumption that the vertical shear on the horizontal velocity fields is negligible, 
which is reasonable for a fluid with shallow hypothesis, (assuming also no rotation for sake 
of simplicity) and using (16. Il l we obtain a 2D momentum equation that reads: 

hi J 

V h dp = -p(w h V hT )(<rdB t ) h , 
V • w = 0, 

V-er = 0, 

V- V-a = 0. 

Here the supperscript, h, indicates an horizontal vector. In order to stick to the no vertical 
shear hypothesis the horizontal velocity components can be seen as averaged along the ver- 
tical axis between an upper free surface, h u , of the fluid and the bottom topography, hb- As 
for the turbulent component it is natural in the same way to rely on a 3D Brownian uncer- 
tainty vector independent of the depth. Its quadratic variation consequently does not depend 
on depth either. This choice ensure a strict respect of a null vertical shear and provides a 
diffusion tensor that is independent from depth. The 3D velocity will be thus defined as a 
velocity field of the form: 

'U{x, y, t)\ fu(x, y, t)dt + {a{x, y, t)d&£) a 
V(x, y, t) = v (x, y, t)dt + (ct(x, y, t)dB?) y | . (6.5) 

W(x, t) ) \ W ( X: t )dt + {a{x, y, t)dB?) 2 



(6.4) 
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Notice that the noise is a 3D white noise defined on the plane and that the vertical velocity 
component depends on the height. Considering the free surface, h u , as a material surface in 
which no fluid crosses the surface, we have, from the stochastic transport principle: 

1 d 2 

d(h u ) t +(Vh u w h -- Y, Q Q {a i jh u ))dt+Biv(h u {trdB t ) h ) = -w{h u )-{adB t ) z , 



(6.6) 



where, w(h u ), denotes the vertical velocity at point (a;, y) of the surface h u and (i,j) h indicates 
summation over horizontal indices. In the same way, for the stationary topographic depth we 
obtain: 

1 d 2 

(Vh b w h - - Y, Q-Q- (aijh b ))dt + Div(h b crdB t )) = -w{h b ) - {crdB t ) z . (6.7) 

The integration of the hydrostatic relations ( 16.1) from a depth z up to the free surface, 
gives for a constant density: 

p'(x,y,z,t) -p'(x,y,h u ,t) = gp(h u - z). (6.8) 

Here enters the two possible choices on the hydrostatic balance. For a hydrostatic balance de- 
fined up to a noise, it is possible to impose an isobaric upper boundary condition to the global 
pressure: p u (t) = p(x, y, h u , t). This condition implies necessarily p'(x, y, h u ,t) = p u (t), 
accompanied with dp(x,y,h u ,t) = as a boundary condition for the random term. This last 
constraint is inappropriate for a global hydrostatic balance as in this case the random com- 
ponent is depth free. In this case, it is necessary to impose a boundary condition only to the 
deterministic pressure function. Let us note that in this latter case, (16.3b imposes to the uncer- 
tainty to be orthogonal to the depth free vertical component gradient d z dp t = w'V T (rdBt = 
0. This also introduces a locally noisy variable pressure on the upper surface which might 
be difficult to control. The former case, is less strict and allows a supplementary degree of 
freedom for the uncertainty. 

In both cases we have nevertheless from i 



p'(x,y,z,t)=gp(h u -z)+p u , (6.9) 

and 

V h p'(x,y,z,t) = g P V h h u (x,y,t), (6.10) 

where we have assumed that p u is horizontally uniform. The momentum equation of the 
horizontal velocity field becomes thus 

f)iii h 1 ^^ r) 2 

(_ + ^V^-- £ (o^)),--^, (6.11) 

and is independent of the vertical coordinate. Integrating now the velocity drift divergence 
free constraint along the z axis , we get: 

- V • w h (h u - h b ) = w(x, y, h u , t) - w(x, y, h bl t). (6.12) 
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In the same way, the integration along depth of the uncertainty divergence free constraint 
provides the relations: 

V ft -cr' l =0, (6.13) 

V h ■ (V' 1 ■ a h ) =0. (6.14) 

The horizontal component of the uncertainty vector is hence divergence free and the horizon- 
tal part of the quadratic variation tensor obeys a 2D uncertainty mass preservation constraint. 
Introducing ( 16.1 3116. 14t in the difference between ( 16. 6i and d6.7t gives for h = h u — hi,: 

dh t + (V • (hw h ) -\Y. d Xt d X] ( ai] h)dt + Vh{(rdB t ) h = 0. (6.15) 

(*,i) h 

The 2D whole shallow water model with uncertainty is hence finally given by 

,dw h 



(-^- + w h V T w h --^ d Xi d Xj { ai jW h )p = -gpVh u , 

(*,i) h 

t + (V • (hw h ) -\Y1 dx&MMdt + Vh{adB t ) h = 0, 



(6.16) 



dh, , 2 

hi 
V h dp = -p(w h V T )((rdB t ) h , 

V ■ cr h = 0, 

V • (V • a h ) = 0. 

The surface is here driven by a stochastic pdes that involves a multiplicative noise. The 
horizontal velocity fields is coupled to this stochastic dynamics through the upper surface 
gradiant. This system is nevertheless unclosed. As a matter of fact, in this inviscid case, the 
second equation cannot be used anymore to defined the diffusion tensor as one can always 
find a random pressure for any given diffusion tensors. The diffusion tensor or its quadratic 
variation must hence be modelled. 

We can point out it is possible to simplify greatly this system if one seeks only a mean 
horizontal velocity field. As a matter of fact, taking the surface conditional expectation, h u , 
upon a given initial condition for the surface provides us the following continuity equation 
for the shallow water model with uncertainty: 

^ + V • (hw h ) - \ Y, d ** d *Mj~ h ) = °> ( 6 - 17 ) 

(hi) h 

This amounts then to solve the system: 

f)iii h 1 

(-gf + w h V T w h - - Y d Xi d Xj ( aij w h ))p = -gpVh u , 

(.i,j) h 

^ + V • (hw h ) \ J2 d Xi d Xj ( ai fh) = 0, (6 - 18) 

(*j) h 
V • (V • a h ) = 0. 

The same system can also be obtained specifying the uncertainties lie on iso-height surfaces 
since in that case the Brownian random terms cancels out. In both cases we get a determin- 
istic formulation of the shallow water equations under uncertainty. The first momentum and 
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the free surface evolution include both a similar subgrid diffusion model accounting for the 
action of the uncertainty term associated to the noisy fluid particle location. Proceeding to 
an integration along the depth direction and writing the system in conservative form we get 
finally: 

(^- + V • {hw h {w h ) T +g-h 2 l) -ij] d Xi d Xj ( aij hw h ))p = -gphVh, 



(*j) h 
dt v ' 2 



' V-{hw h )-\ ^^^(a^) = 0, (6 ' 19) 



V • (V • a h ) = 0. 

A stochastic system of the same conservative form could as well be obtained from the spde's 
( 16.16l >. The derivation applied here is general, and could be extended to the different forms 
of geophysical flow models. 

Another example of the eventual benefits brought by the application of our formalism 
concerns the establishment of reduced order dynamical systems. We briefly describe such a 
derivation in the following section. 

7. Application to low order dynamical system modeling. The constitution of low or- 
der dynamical models to describe the evolution of fluid flows arouses an intense research in- 
terest in several domains ranging from climatic study [|25l . flow control ]3]|29] or atmospheric 
sciences 1371 [131 . Reduced order model are usually defined from a Galerkin projection on a 
reduced basis specified from experimental measurements or numerical simulation data. One 
of the simplest way to define such empirical basis functions stems from the Karhunen-Loeve 
decomposition (referred as proper orthogonal decomposition - POD - in the fluid mechanics 
domain). In this decomposition the basis functions encodes the direction of higher variance 
and are solutions of eigenvalues problems associated either to the two points correlation ten- 
sor or to the temporal correlation tensor (see for instance [fTSl ). This expansion is easy to 
implement, and leads through a Galerkin projection to a system of ordinary differential equa- 
tions that enables representing with only few modes the evolution of complex flows. Never- 
theless, the system built reveals often unstable in practice and does not allow operating long 
term forecasts even if it fits perfectly to the data JT]|30]. To mitigate structural instabilities, 
artificial damping terms are usually introduced in the reduced model. One of the classical 
solution consists in adding a constant viscosity acting in the same way on all the POD modes 
l35l . Approaches involving different viscosity values for the set of modes have been also 
considered 1341 and spectral vanishing dissipative model |fl6l have been introduced. Those 
damping terms added to the flow kinematic viscosity enables improving the system's numer- 
ical stability but require often a proper tuning of the set of parameters involved. Some of the 
methods proposed may even require fixing a priori a great number of parameters. The prin- 
ciple followed by all those techniques consists in artificially reintroducing a lost dissipative 
mechanism attached to the mode discarded by the severe truncation operated at the Galerkin 
projection stage. The inclusion of uncertainty in the flow dynamics brings naturally such a 
mechanism, with the additional advantage that the parameters can be properly fixed from the 
data. In the following section we recall briefly the principles of the proper orthogonal de- 
composition, and we show how a Galerkin projection of the Navier-Stokes equations under 
uncertainty onto a reduced set of empirical modes enables to get a dynamical system of low 
order. 

7.1. POD basis. POD decomposition has been widely used by different authors as a 
technique to obtain approximate descriptions of the large scale coherent structures in laminar 
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and turbulent flows. Given an ensemble u(x, ij) of velocity fields obtained experimentally at 
M different discrete instants, POD provides a set of M mutually orthogonal basis functions, 
or modes, <j>i{x), which are optimal with respect to the kinetic energy. 

Considering such a decomposition enables to write the velocity field as an average u 
with fluctuations captured by a finite set of modes: 

M 

u(x,t) = u + J2 b i( t )Mx), (7-1) 



where, /, denotes a temporal average of function /. Denoting by (, ) the inner product asso- 
ciated to the finite energy functional space L 2 (fl) and |/|n = (/, f) 1 ^ 2 its induced norm: 

(u,ip) n = / uij>ds, (7.2) 

seeking an optimal representation subspace for u(x, U) on the domain il and along the sam- 
pling time comes to find an ensemble of functions that maximize: 



\(u,ip)\ 2 2 , subject to {i/),i/>)a = 1, (7.3) 

The solution, <fi, of this constrained optimization problem satisfies the following eigenvalue 
problem: 

/ K(x,x')(j> k (x)dx = (j>k(x)\ k , (7.4) 

Jn 



with 



1 *r 

K(x,x') = (u(x,t)u(x',t)) = — y^ u(x, tj)u(x', U) 



i = l 



The correlation tensor is a linear, positive semi-definite and self-adjoint. As a consequence, it 
admits by the Hilbert-Schmidt theorem a spectral representation, which guaranties a solution 
to problem (I7.3b . The eigenvalues are real and positive and by the Mercer theorem, the spatial 
autocorrelation can be represented as the (uniformly convergent) series: 



+oo 



R(x, x') = J2 h<t>k{x)4>l{x f ). (7.5) 

i=0 

The spatio temporal symmetry of the representation in terms respectively of the temporal 
{hi, i = l,..., M} or spatial {(pi, i = 1, . . . , M} expansions allows interchanging the two 
point correlation tensor with the two times correlation tensor(J38). Such a procedure is numer- 
ically advantageous when the temporal dimension is lower than the spatial dimension. This 
new eigenvalue problem, whose eigenvectors are the temporal coefficients, bk, is obtained by 
expressing in (17. 4t the spatial modes as a linear combination of the velocity fields: 

1 



<t>k{x) = —b k (t)u(x,t) 
Afe 

The eigenvalues of this new problem are identical to those obtained for the spatial modes. 
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7.2. Galerkin projection. A Galerkin projection enables to rewrite a partial differen- 
tial equation (PDE) system as a system of ordinary differential equation (ODE). According 
to this procedure, the functions which define the original equation are projected on a finite 
dimensions subspace of the phase space (in this case, the subspace generated by the first s 
modes). 

A Reynolds decomposition of the Navier-Stokes equation under uncertainty can be for- 
mulated by separating the deterministic drift component into a mean w and a fluctuating w' 
parts: w'(x,t) = w(x, t) — w(x). It describes the evolution of the fluctuating deterministic 
part of the velocity fields. This system reads: 

vy 

d t w' + w'Vw + wVw' + w'Vw' — w'Vw' -\ i/A(w + w') — 



- y^d Xi d x . (aij(w' + w)) = 0, 



2 

hi 

Vdpt = -p(wV T )crdBt + /j,AcrdB t , 

V-™ = 0,V-cr = 0,V-(V-a) = 0. (7.6) 

Let us assume now that the deterministic drift component, w, lives in a subspace of finite 
dimension spanned by a set of basis functions $ = {</>i,i = 1, . . . , m} determined from a 
sequence of observed motion fields {u(- ■ ■ ,tj),j = 1, ...,/} much longuer than the number 
of empirical modes (m < /). The highly random oscillating terms crdB t are assumed to live 
in the complement orthogonal space spanned by the basis $ cs = {4>i, i = m + 1, . . . , /}. 
Applying a Galerkin projection of system ( 17.6) onto the truncated basis, <1>, we get 

Vp' 

dtw' + w'Vw + wVw' + w'Vw' — w'Vw' -\ i/A(w + w') — 



\ jy* t d Xj (cnj K + w)), $i I = 0. (7.7) 



Rewriting d7.7| i in terms of the POD temporal coefficients, we obtain a quadratic system of 
ODEs. For every j < m modes, the system reads: 

77 m s m 

—jT = ik + ^2kkh + *^2^b l Ci j kbj Vfc = 1, ...,m (7.8) 

i— 1 i— 1 j—i 

where 



wW4>i4>jdx+ I (pi'Vwcpjdx- 
n 



(7.9) 



o / y^Pxk d xe( a kt<Pi)<l>jdx - / —Acjjicjjjdx, 

"ill,// "Si 

(kjk = / 4>jV(t>i(t> k dx, (7.10) 

Jn 

h = / Vp'4> k dx - — / Aw4> k dx- 
Jn Re J n 

- / Y]dxj d X;l (aijw)(l> k dx - y~] Xj / d>jV<pj(l> k dx 
1 Jn „■ " „_, Jn 



(7.11) 



/Q i,3 3=1 
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All those expressions capture different physical phenomenon. The first linear term (17.9) 
describes the interaction between the mean flow and the fluctuating field. It also includes 
viscous effects of the resolved modes and dissipation caused by the uncertainty component. 
Nonlinear advection effects are reported by (17.10) . The independent term ( 17.11) takes into 
account mean flow dissipation due both to friction and to the action of the random uncertainty, 
convective contribution of the modes and the pressure field influence. Boundary conditions 
and symmetry make the pressure term vanish in the particular case of wake flows. As a matter 
of fact, each mode satisfies the continuity equation, to give: 

/ Vp'<t> k dx = f p'tpkds, 
Ja Jc 

where C is the boundary of domain ft. Works of (8) and |28l demonstrated that for wake 
flow configuration, the latter expression is negligible compared to the other terms. Direct 
calculation of each term of the system ( 17.8) can be avoided by relying on polynomial iden- 
tification or optimal control techniques [1]. The later provide in addition a way to estimate 
an adapted initial condition of the system and to consider eventually an error term on the 
dynamics. Least squares identification or variational assimilation techniques HI , hides the 
necessity of an additional disipation mechanisms as it is assumed to be carried by the data. 
However, even in that case it is usefull to formulate in an explicit way the precise form of the 
uncertainty involved. 

The simplest procedure to specify the quadratic variation tensor, a, associated to the 
uncertainty terms, consists in identifying it to the measurements variance: 

m \ / m 

u(x,U) - u -y^bj(tj)ct>j(x) u(x,U) - u -y^bj(tj)4>j(x) 

3=1 J V 3 = 1 

Here it is assumed that the measurement variance tensor respects the incompressibility con- 
straint V • V • a = 0. A strict respect of this constraint would require an additional projection 
defined through a constrained least squares fitting. However it can be noticed, that the sim- 
ple scheme above provides a very intuitive scheme as the diffusion process of the large scale 
flow component is proportional to the residuals variance tensor. Locally the greater the dis- 
crepancy between the reduced model and the data the stronger the diffusion of the velocity 
components. 

This procedure provides immediately an expression of the covariance tensor associated 
to the uncertainty in terms of the complement space basis functions: 

/ 
E(a(x)dB t <j(x')dB? ) = J2 ^MxWk(x'). (7.12) 

k— m+1 

The uncertainty term can then be classically written in a spectral form as: 

/ 
cr(x)dB t = J2 ^) /2 <t>dP](t), (7.13) 

j— m+1 

where f3j(t) denotes a set of uncorrelated scalar standard Brownian motions. This represen- 
tation supplies thus a simple scheme for the uncertainty sampling and hence allows random 
drawings of the whole flow field: 

m f 

?'— 1 j— m+1 




22 E. Memin 

Such a formulation should enables to set up stochastic filtering procedures to estimate state 
variable of interest from partial observation like satellite images J4| in which the reduced 
order modeling would have the great advantage to offer the possibility of an immediate sig- 
nificant augmentation of the ensemble members. 

8. Conclusion. In this paper we have proposed a decomposition that allows us mod- 
elling the action of uncertainties related to the fluid particles displacement on a resolved 
component of the flow dynamics. The random uncertainty fields encodes numerical artefacts 
or unresolved physical processes that have been neglected in the momentum balance. They 
are defined through diffusion tensors that have to be estimated or specified. In the former 
case, the estimation can be led from a set of vectorial Poisson equations and two divergence 
free constraints. In the latter case the system reduces to a Navier-Stokes equation with a sub- 
grid modelling defined from an anisotropic diffusion that involves only the variance of the 
random term together with a global constraint on this variance. This formulation gives a way 
to define appropriate subgrid diffusions from the flow uncertainties specification and might 
be useful to build large scale stochastic climatic or geophysical models. It brings also a new 
point of view to some usual subgrid modelling and introduce additional constraints on those 
models. The relation obtained between the anisotropic diffusion and uncertainties should be 
particularly useful for ensemble data assimilation methods and ensemble forecasting. 
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